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ABSTRACT 


J 

In  Part  I,  reduced  basis  techniques  are  applied  to  the  prol-  -•«  of  the  nonli¬ 
near  analysis  of  the  dynamics  of  unrestrained  flexible  structures.  The  reduced 
bases  used  consisted  of  mode  shapes  of  the  structure  as  well  as  some  modal  deri¬ 
vatives.  The  technique  was  tested  on  a  simple  spacecraft  structure.  The 
numerical  results  indicated  that  the  technique  did  not  appear  very  promising  for 
this  type  of  problem. 

In  Part  II,  a  finite  element  technique  is  used  for  analysis  of  very  flexible 
structures  undergoing  deployment  maneuvers.  The  structure  is  assumed  to  consist 
of  flexible  bars  attached  to  a  rigid  mass.  The  description  of  elastic  defor¬ 
mation  is  based  on  the  total  Lagrangian  formulation  which  allows  finite  rota¬ 
tion.  Numerical  tests  demonstrates  the  validty  of  the  present  approach. 


A  REDUCED  BASIS  TECHNIQUE  FOR  THE  NONLINEAR  ANALYSIS 
OF  THE  DYNAMICS  OF  UNRESTRAINED  FLEXIBLE  STRUCTURES 


1.  INTRODUCTION 

For  large  lightweight  spacecraft  which  are  undergoing  certain  types  of 
manuevers  such  as  slewing  motions  for  retargeting  or  rapid  large  angle  reorien¬ 
tations,  large  elastic  deformations  as  well  as  large  rigid  body  motions  can 
occur.  These  large  deformations  will  be  coupled  with  the  rigid  body  motions 
resulting  in  a  complicated  set  of  nonlinear  differential  equations  of  motion. 
Applying  a  finite  element  approximation  to  such  a  problem  for  a  complicated 
structure  results  in  a  system  of  equations  involving  a  large  number  of  degrees 
of  freedom  [l].  In  order  to  avoid  using  excessive  amounts  of  computer  time,  it 
is  desirable  to  reduce  the  number  of  unknowns  using  some  sort  of  reduction 
method.  One  method  is  to  apply  the  classical  Rayl ei gh-Ritz  approximation  to  the 
system  of  matrix  equations  resulting  from  the  finite  element  discretization  [ 2j 
thus  combining  the  modeling  versatility  of  the  finite  element  method  with  the 
reduction  in  variables  due  to  the  Raylei gh-Ritz  method.  The  accuracy  of  the 
reduction  method  depends  of  course  upon  the  choice  of  Rayleigh-Ritz  approxima¬ 
tion  functions  or  basis  vectors.  According  to  [3],  these  basis  vectors  should 
satisfy  the  following  criteria: 

(1)  The  vectors  must  be  linearly  independent  and  span  the  space  of  solu¬ 
tions  in  the  neighborhood  of  the  point  considered  on  the  solution  path. 

(2)  Their  generation  should  be  simple  and  computationally  inexpensive. 

(3)  They  must  have  good  approximation  properties  in  the  sense  that  they 
provide  highly  accurate  solutions  on  a  large  interval  of  the  solution 


For  the  solution  of  static  problems,  Noor  [2,3]  has  recommended  the  use  of  a 
nonlinear  solution  and  its  various  order  path  derivatives  as  basis  vectors. 

Chan  and  Hsiao  [4]  have  investigated  the  use  of  a  normalized  set  of  correction 
vectors  derived  during  the  course  of  a  nonlinear  analysis  step  as  the  basis  vec¬ 
tors  for  geometrically  as  well  as  materially  nonlinear  static  problems.  For 
nonlinear  dynamic  problems,  modal  superposition  techniques  are  a  popular  reduc¬ 
tion  method  [5-11].  These  techniques  are  based  in  part  upon  the  principle  of 
local  mode  superposition  [12]  which  states  that  small  harmonic  motions  may  be 
superimposed  upon  large  displacements  and  that  small  forced  motion  may  be  repre¬ 
sented  in  terms  of  the  nonlinear  (tangent  stiffness)  frequency  spectrum. 

Nickel  1  [5]  and  Resmeth  [6]  have  used  a  modal  reduction  technique  in  which  the 
modes  and  frequencies  are  updated  at  each  time  step  in  the  solution.  Tnis 
approach  requires  a  repeated  solution  of  the  free  vibration  eigenvalue  problem. 
Other  investigators  [7-11]  suggest  using  the  initial  eigenvalues  throughout  the 
analysis.  The  nonlinearities  must  then  be  taken  into  account  by  a  psuedo  force 
approach  in  which  the  nonlinear  terms  are  included  in  the  force  terms  on  the 
right  hand  side  of  the  equation.  Bathe  and  Graceuski  [11]  recommend  this 
approach  for  problems  with  small  or  local  nonlinearities.  Noor  [3]  has  reported 
favorable  results  using  a  combination  of  mode  shapes  expanded  about  the  unloaded 
equilibrium  configuration  and  a  set  of  mode  shapes  expanded  about  the  nonlinear 
static  solution.  Wilson  et.  al  [13,14]  have  used  a  sequence  of  orthodgonal  Ritz 
Vectors  as  a  reduced  basis  for  linear  dynamics  problems.  Their  approach 
requires  that  the  loading  terms  be  factored  into  a  sum  of  space  vector 
multiplied  by  time  functions. 

Idelsohn  and  Cardona  [15,16]  have  shown  that  for  highly  nonlinear  dynamic 
systems  such  as  vibrating  cantilevered  beams  that  modal  superposition  techniques 
require  an  updating  of  the  basis  too  frequently  and  lead  to  inaccurate  results. 
They  reported  improved  results  with  not  only  the  tangent  eigenmodes  but  also 
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some  derivatives  of  them  with  respect  to  the  modal  amplitude  parameters  as  the 
basis  vectors. 

All  of  the  work  mentioned  above  has  been  involved  with  fully  restrained 
structures.  The  objective  of  this  report,  therefore,  is  to  study  the  use  of 
reduced  basis  techniques  in  the  solution  of  problems  dealing  with  nonlinear 
unrestrained  flexible  structures  such  as  those  studied  in  [lj. 


2.  MATHEMATICAL  FORMULATION 

2.1  Equations  of  Motion 

The  motion  of  an  unrestrained  flexible  body  can  be  described  by  the 
following  three  sets  of  equations: 

1.  The  conservation  of  linear  momentum 

2.  The  conservation  of  angular  momentum 

3.  An  equation  describing  the  deformations 

The  finite  element  equations  of  motion  for  an  unrestrained  deformable  body 
undergoing  large  deformations  have  been  derived  in  [lj.  The  equations  evaluated 
at  time  tn  are 

uT  +  £E(3n)  -  F(fl,a>  (2.1) 

The  modal  displacements  including  the  rigid  body  ones  are  denoted  by  3.  The 
mass  matrix  is  M  and  F  is  the  applied  force  term.  The  vector  of  internal  modal 
forces  is  F^g)  which  is  a  nonlinear  function  of  the  modal  displacements.  Eq. 
(2.1)  is  nonlinear  and  can  be  solved  by  a  modified  Newton-Raphson  interations. 
Letting  subscript  i  represent  the  ith  iteration,  eq.  (2.1)  can  be  written  as 


n*.n  „  .  n  _n 

5!  S1+1  +  &rAfli  =  Et-j 


where 


.  n  n  n 

ASi  =  3i+i  -  9i 


(2.2) 


(2.3) 


’  A 


4 


(2.4) 


r-  n  — .  n  _  n 

£t(  ■  Ei  -  Es, 

Kj  is  the  tangent  stiffness  matrix  and  F$  is  the  initial  stress  force  vector. 

Eq.  (2.2)  can  be  integrated  timewise  using  an  implicit  scheme  based  on  the  tra- 
pizoidal  rule  [l].  For  a  large  system,  the  cost  of  this  computation  can  be 
substantial.  Reduced  basis  techniques  will  be  applied  to  Eq.  (2.2)  in  an  effort 
to  reduce  those  costs. 

2.2  Reduced  Equations 

The  Rayleigh-Ritz  technique  is  applied  to  Eq.  (2.2)  by  approximating  3"  as 
follows: 


(2.5) 


where 


£  =  [jfei  *2  •••  A™] 


(2.6) 


r  is  a  matrix  made  up  of  m  linearly  independent  basis  vectors  <j^ ,  •••»  &111 

where  m  is  much  less  than  the  total  number  of  degrees  of  freedom  for  the 
problem.  The  vector  Xn  represents  m  undetermined  coefficients  which  must  be 
solved  for  from  the  rduced  equations.  From  Eq.  (2.5), 


*n  ;n 
Hi  -  £ 


(2.7) 


Ci  =  £  Ci 


(2.8) 


Substituting  Eqs.  (2.5),  (2.7)  and  (2.8)  into  Eq.  (2.2)  and  multiplying  through 
by  rT  results  in 


no  "n  -tTi  un  — n 

n  Xi+i  +  =  ft  i 


(2.9) 


where 


M*1 


rT  Mn  r 


(2.10) 


K " 

~  I 


r 


(2.11) 


T?.  =  rT  Fj . 

1  ]  ~  /V  I  ] 


(2.12) 


Eq.  (2.9)  is  the  reduced  set  of  m  equations  which  are  solved  timewise  by  the 

same  scheme  as  in  [lj.  Note,  however,  that  since  FT  =  jFy(jg ,jg ) ,  it  is  necessary 

to  transform  back  to  the  full  coordinates  via  Eq.  (2.5)  after  each  iteration  in 

order  to  recalculate  FT  for  the  next  iteration. 

~  l 

3.  SELECTION  OF  REDUCED  BASIS  VECTORS 

The  most  difficult  and  critical  part  of  the  reduced  basis  technique  is 
choosing  what  vectors  to  use  as  a  basis.  Two  different  choices  of  vectors  were 
considered  in  this  study.  One  choice  consisted  of  a  modal  method  in  which  the 
basis  vectors  consisted  of  the  system  eigenvectors.  For  the  other  choice,  some 
eigenvector  derivatives  were  added  to  the  basis.  Both  of  these  methods  are 
described  in  more  detail  below. 


3.1  Modal  Method 

The  use  of  this  method  is  very  similar  to  the  use  of  modal  superposition  for 
a  linear  problem.  The  basis  vectors  consist  of  the  first  few  modes, 

r  =  $  (3.1) 

where  $  is  an  Nxm  matrix  consisting  of  the  m  eigenvectors  ...,  £m,  that 


V  VS  V"  X.T  T^rr^rr*  v  ^tt*V  rJr?'rZv~Zli~.  *'“  :^r>" \  ’J". 

Unlike  the  linear  problem,  however,  the  Ky  and  M  matrices  will  not  in  general  be 
diagonal  since  the  Ky  and  My.  matrices  are  changing  with  time.  The  eigenvectors 
are  found  by  solving  the  linearized  eigenvalue  problem. 


(Ky  -  u>2Mn)  tr  =  o  r=l,2,...m  (3.3) 

The  eigenvectors  can  be  found  from  Eq.  (3.3)  using  any  method.  Since  the 
problems  in  this  study  didn't  involve  a  large  number  of  degrees  of  freedom,  the 
Jacobi  method  was  used  to  calculate  all  the  modes.  For  larger  problems,  more 
efficient  schemes  such  as  subspace  iteration  or  Lanczos  method  should  be  used. 


3.2  Modal  Derivatives  Method 

Because  the  problem  is  nonlinear,  the  actual  eigenvalues  and  eigenvectors 
will  be  functions  of  the  displacements  or  alternatively,  they  will  be  func¬ 
tions  of  the  modal  participation  factors  z.  Thus, 

Ar  ■  (3*4) 

and 

m 

A%  =  l.  trWlr  (3-5) 

r=l 

If  Eq.  (3.4)  is  expanded  in  a  Taylor  series  around  the  reference  position  of  c)^. 


m  3(fi 


’r  =  Ar  +  Azs  + 


s-1  >9zs 


(3.5) 


In  the  modal  method  it  is  the  first  term  of  the  Taylor  series  which  is  used  as  a 
basis.  In  the  modal  derivative  method  some  of  the  modal  derivatives  are  added 
to  the  basis  so  that 


I  =  [Ai  •••  A, 


3<J> 


1 


9$ 


m  **•  3z$ 


(3.6) 


We  needn't  add  all  of  the  modal  derivatives  to  the  basis,  the  choice  of  which 


ones  to  keep  will  depend  upon  the  particular  problem  being  solved.  By  including 
these  extra  terms  from  the  Taylor  series  in  the  basis  a  more  accurate  represen¬ 
tation  of  the  response  may  be  obtained. 

The  modal  derivatives  can  be  found  by  differentiating  Eq.  (3.3)  with  respect 


to  the  modal  amplitude  z$  [16]  and  then  solving  the  following  set  of  equations 
[17]: 


where 


T  3Kt  3Kt 

~rs  -  ^r^r  3z  ~r^  ~  Tz7  ~r 


(3.8) 


The  *  denotes  the  matrix  and  vector  obtained  by  deleting  the  kth  row  and  column 


where  k  is  selected  such  that  the  absolute  value  of  the  kth  component  of  is  a 
maximum.  The  derivative  of  the  stiffness  matix  is  found  as  follows  [ 16]  : 


n  o  r\ 

l  4>s  i  a ~ 
j-1  ~  J 


(3.9) 


Thus  in  order  to  determine  the  modal  derivatives  it  is  necessary  to  calculate 
the  derivative  of  the  stiffness  matrix  with  respect  to  all  the  modal  displace¬ 
ments. 


4.  NUMERICAL  RESULTS 

Both  the  modal  method  and  the  modal  derivative  method  is  applied  to  the 
sample  problem  illustrated  in  Fig.  1.  A  simple  spacecraft  consisting  of  a  rigid 
mass  and  a  flexible  beam  appendage  is  subjected  to  a  sinusoidally  varying  moment 
applied  to  the  rigid  mass.  The  frequency  of  the  applied  moment  is  equal  to  the 
first  bending  frequency  for  the  undeformed  structure.  The  beam  was  modelled  by 
two  of  the  three  node  beam  elements  developed  in  [lj  which  resulted  in  13  total 
unknown  degrees  of  freedom,  including  the  rigid  body  rotation  angle  of  the  rigid 


mass.  The  equations  of  motion  were  integrated  using  a  time  step  size  of  0.01 
sec.  The  numerical  results  for  each  case  are  presented  below. 

4.1  Modal  Method 

The  results  using  various  number  of  modes  in  the  reduced  basis  are  shown  in 
Fig.  2  where  the  lateral  tip  displacement  of  the  beam  is  plotted  as  a  function 
of  time.  The  first  case  considered  uses  five  modes  in  the  basis:  One  rigid 
body',  the  first  two  bending  modes  and  the  first  two  axial  modes.  In  the  second 
case  considered  the  basis  consisted  of  nine  vectors:  the  rigid  body  mode  plus 
four  bending  and  four  axial  modes.  The  original  basis  was  used  throughout  the 
computations  for  these  two  cases.  The  third  case  considered  used  three  modes, 
one  rigid  body  and  two  bending  modes  as  the  reduced  basis.  For  this  case, 
however,  the  three  modes  were  recomputed  at  each  time  step. 

From  Fig.  2  it  can  be  seen  that  in  none  of  the  non-updated  cases  were  the 
reduced  equations  sufficient  to  follow  the  true  response  of  the  complete  system. 
The  more  modes  included  in  the  basis,  the  better  the  results,  but  eventually  the 
results  become  inaccurate,  even  with  nine  modes  used  in  the  basis.  The  results 
were  much  better  when  the  modes  were  updated  at  each  time  step,  but  this 
entailed  a  great  deal  of  extra  computational  expense.  Also,  the  truncation 
errors  due  to  repeated  change  of  basis  accumulated  so  that  eventually  the 
results  even  for  this  case  were  in  error. 

4.2  Modal  Derivative  Method 

In  addition  to  using  some  modes  in  the  reduced  basis,  some  modal  derivatives 
are  also  added.  The  results  for  these  cases  are  pictured  in  Fig.  3.  The  first 
case  considered  included  a  total  of  five  basis  vectors  -  one  rigid  body  mode, 
the  first  two  bending  modes  and  the  derivatives  of  the  first  bending  mode  with 
respect  to  its  own  amplitude  and  with  respect  to  the  second  bonding  modes  ampli¬ 
tude.  The  second  case  considered  included  9  basis  vectors  -  the  first  five 


modes  (one  rigid  body  +  first  four  bending)  plus  four  modal  derivatives,  the 
derivatives  of  the  first  bending  mode  with  respect  to  the  amplitudes  of  bending 
modes  one,  two,  and  three  plus  the  derivative  of  the  second  bending  mode  with 
respect  to  its  own  amplitude. 

From  Fig.  3,  however,  it  does  not  appear  that  the  addition  of  the  modal 
derivatives  has  improved  the  results  very  much.  In  fact,  by  comparing  Fig.'s  2 
and  3,  it  appears  that  the  results  depend  upon  the  total  number  of  basis  vectors 
included  rather  than  upon  their  nature. 

5.  CONCLUSIONS 

From  the  numerical  results  obtained  in  this  study  it  appears  that  the  use  of 
a  reduced  basis  technique  for  the  solution  of  the  type  of  problems  studied  here 
is  not  very  promising.  It  may  be  possible  that  a  different  choice  of  basis  vec¬ 
tors  would  give  improved  results. 

The  best  results  were  obtained  by  updating  the  modes  at  each  time  step. 
However,  this  procedure  involves  considerable  computational  expense.  In  addi¬ 
tion,  the  truncation  error  buildup  due  to  the  repeated  change  of  basis  even¬ 
tually  led  to  large  errors. 

It  appears  for  this  type  of  structure  modelled  by  the  type  of  elements  used 
here  that  the  changes  in  the  modal  properties  are  large  enough  that  reduced 
basis  techniques  are  not  practical  due  to  the  large  number  of  basis  updatings 
necessary  to  accurate  model  the  response. 
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PLUS  MODAL  DERIVATIVES  INCLUDED  IN  BASIS 


FINITE  ELEMENT  MODELING  OF  DEPLOYMENT 
DYNAMICS  OF  SPACE  STRUCTURES 

1.  INTRODUCTION 

The  objective  of  the  research  project  has  been  to  develop  an  efficient 
finite  element  analysis  technique  for  study  of  deployment  and  retraction  dyna¬ 
mics  of  large  flexible  space  structures  in  earth  orbit. 

With  the  advent  of  the  space  shuttle,  it  has  become  feasible  to  construct  a 
large  space  structure  which  can  be  used  as  permanent  communication  and  obser¬ 
vation  posts.  However,  due  to  its  enormous  size,  a  typical  space  structure  can¬ 
not  fit  into  the  cargo  bay  of  the  current-generation  space  shuttle  in  fully 
deployed  mode.  Instead  it  must  be  folded  into  a  much  smaller  volume  for  effi¬ 
cient  stowage.  Once  placed  in  earth  orbit,  the  space  structure  is  separated 
from  the  shuttle  and  deployed  into  full  dimension  for  operational  use.  Also  full 
or  part  of  the  space  structure  may  undergo  a  series  of  deployment  and  retraction 
manuevers  as  a  part  of  normal  operation.  Large  change  in  r-ometrica'1  con¬ 
figuration  of  space  structures  during  deployment,  retraction  and  other  maneuvers 
implies  large  change  in  inertia  properties  and  coupling  between  gross  rigid  body 
motion  and  elastic  deformation.  In  particular  space  structures  are  extremely 
flexible  due  to  minimum  weight  requirement  and  thus  are  likely  to  experience 
large  elastic  deflection. 

Design  of  space  structures  requires  a  very  accurate  and  reliable  analysis 
technique.  This  is  due  to  the  fact  that  extremely  large  size  of  space  struc¬ 
tures  and  zero  gravity  environment  make  any  full-scale  test  on  ground  prac¬ 
tically  impossible.  Therefore  more  reliance  on  small-scale  testing  and 
theoretical  analysis  is  inevitable.  Due  to  its  versatility  in  modeling  compli¬ 
cated  geometry,  finite  element  method  is  the  most  widely  used  tool  for  analysis 
of  various  structural  dynamics  problems.  Especially,  for  such  a  highly  nonli- 


near  time-dependent  problem  as  deployment/retraction  dynamics  of  space  struc¬ 
tures  with  large  change  in  inertia  and  stiffness  properties,  finite  element 
method  seems  to  be  the  only  viable  tool  for  theoretical  analysis. 

There  have  been  many  papers  discussing  the  methods  for  construct!  iy 
equations  of  motion  governing  the  deployment  of  a  various  types  of  flexible 
appendages  from  a  spacecraft  [1-9].  The  main  objectives  were  the  investigations 
of  the  effects  of  deploying  flexible  appendages  on  the  attitude  motion  and  sta¬ 
bility  of  a  spacecraft.  However,  in  those  papers,  the  deployment  means  the  only 
extension  of  flexible  appendages  from  the  satellite  in  earth  orbit  and  the 
nonlinear  effects  of  elastic  deformations  were  not  considered.  For  the 
deployment  of  packaged  structures,  Housner  [10]  presented  a  finite  element  for¬ 
mulation  for  transient  analysis  of  two  dimensional  beam-frame  configuration  in 
space. ‘  In  his  paper  a  convective  coordinate  system  was  used  for  each  element 
because  the  Euler  beam  theory  used  in  the  formulation  cannot  handle  large 
elastic  deflection  unless  an  updated  coordinate  system  is  used  for  individual 
elements.  Moreover,  Housner's  formulation  does  not  seem  to  be  convenient  for 
deployment  in  three  dimensional  space. 

In  this  report,  a  finite  element  technique  which  can  handle  deployment  of 
flexible  structures  in  three  dimensional  space  is  presented.  The  structures  is 
assumed  to  consist  of  a  flexible  bars  attached  to  a  rigid  mass.  The  elastic 
deflection  of  bars  can  be  very  large  and  is  coupled  with  gross  rigid  transla¬ 
tional  and  rotational  motions.  The  description  of  elastic  deformation  is  based 
on  the  total  Lagrangian  formulation.  Large  deflection  or  rotation  with  respect 
to  a  body  axis  is  handled  by  introducing  rotational  angles  as  independent 
variables  for  rotation  of  cross-sections.  The  equations  of  motion  are  solved 
timewise  by  a  direct  numerical  time  integration  scheme  using  the  trapezoidal 
rule  and  the  Newton  Raphson  iteration.  The  parametric  studies  have  been  per¬ 
formed  to  investigate  the  effect  of  deployment  on  the  rigid  body  rotational 


motion.  The  effect  of  gravity  gradient  is  also  considered  for  the  deployment  of 
non-symmetric  space  structure  in  earth  orbit. 

2.  EQUATIONS  OF  MOTION 

In  order  to  predict  the  motion  of  spacecraft  with  flexible  bodies  in  earth 
orbit,  various  methods  have  been  developed.  In  such  a  motion,  large  elastic 
deformations  are  coupled  with  the  rigid  body  motions.  The  principles  of  conser¬ 
vation  of  linear  and  angular  momentum  are  adopted  to  describe  the  gross  transla¬ 
tional  and  rotational  motion.  Also  to  describe  the  elastic  deformation  the 
principle  of  virtual  work  is  used.  The  resulting  equations  of  motion  for 
elastic  structures  undergoing  gross  translational  and  rotational  motion  as  well 
as  small  or  large  elastic  deformation  become  a  complicated  set  of  nonlinear  dif¬ 
ferential  equations.  Introducing  finite  element  approximation  the  equations  of 
motion  can  be  written  symbolically  as  follows; 

M  a  +  =  £(a»a)  (2-D 

Note  that  g  denotes  the  nodal  displacements  including  the  gross  translational 
and  rotational  displacements  of  body  axes  and  elastic  displacement  with  respect 
to  the  body  axes.  The  dot  (•)  means  the  time  derivative.  The  details  of  M,  F0, 
and  F  matrices  are  derived  in  reference  15. 

Equation  (2.1)  may  represent  any  type  of  structures.  However,  in  the  pre¬ 
sent  study,  we  are  interested  in  the  development  of  space  structures  which  con¬ 
sist  of  straight  beams  connected  by  hinge  joints.  The  beam  is  assumed  to 
experience  small  strain  but  large  rotations  are  allowed.  This  is  done  by 
allowing  transverse  shear  deformation  and  adopting  the  total  Lagrangian  descrip¬ 
tion  for  the  geometry  and  kinematics  of  deformation.  The  cross-section  of  the 
beam  is  circular  therefore  the  cross-section  is  free  of  warping.  As  stated  in 


the  previous  section,  with  the  present  formulation,  the  kinematics  of  beam 
deformation  can  be  described  with  respect  to  a  single  body  axis  system.  The 
particular  element  used  in  the  present  study  is  the  three  node  element  as  shown 
in  Fig.  1.  Each  node  has  three  translational  displacement  and  three  rotational 
angles  as  nodal  degrees  of  freedom.  Detailed  formulation  is  given  in  Reference 
15. 

3.  NUMERICAL  INTEGRATION  SCHEME 

3.1  Time  Integration  Method 

There  are  many  methods  which  can  be  used  for  numerical  time  integration  of 
Eq.  (2.1)  [13,14].  Generally  these  methods  are  classified  as  either  explicit  or 
implicit  schemes.  Explicit  schemes  are  computationally  simple  but  the  time  step 
size  is  limited  by  stability  consideration.  Implicit  schemes  require  more  com¬ 
putation  per  time  step  but  time  step  size  limitation  are  much  less  stringent. 

In  the  present  study  an  implicit  scheme  based  on  the  trapezoidal  rule  has  been 
used. 

Assuming  all  the  variables  in  Eq.  (2.1)  are  known  at  time  t  =  tn,  the 


equation  at  time  tn+j  = 

3ntl 

tn  +  At  is  written  as 
••n+1  ,  cn+l  _  cn+l 

3  +  £e  "  Z 

(3.1) 

where 
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Using  trapezoidal  rule 


•n+1  *n  ,  At  /**n  ,  ••n+lx 
S  *  3  +  T  (9  +  a  ) 


(3.7) 


Combining  Eq.  (3.6)  and  Eq.  (3.7)  results  in 
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Substituting  Eq.  (3.8)  into  Eq.  (3.1)  yields 


M"+1  [  «  (antl  -  a")  -  4  a"  -  ?]  *  z";1  ■  f1*1  (3.9) 


(At) 


Eq.  (3.9)  is  a  nonlinear  equation  and  can  be  solved  by  the  Newton-Raphson  itera- 

A.  L- 

tion.  Letting  subscript  i  denote  thr.  itn  iteration  number,  Eq.  (3.9)  becomes  as 
follows; 


un+l  r  4  ,_n+l  „n\  4  *n  **ni  x  r  /_n+l\  _  rn+l 

~i  <3(91  '  3  >  ‘It  3  -  3  1  +  £e<3i+1)  =  E,  (3-10> 


As  given  in  reference  15,  F  can  be  decomposed  into 

Ee  *  Ef;1  ‘s"*1  *  Es*1 


(3.11) 


where 


.qn+1  _  n+1  n+1 

A~i  -  a1+i  -  3i 


Note  that  the  Kj  is  the  tangential  stiffness  matrix  and  Fs  is  the  initial  stress 
force  vector.  Using  Eq.  (3.11),  Eq.  (3.10)  can  be  written  as  follows; 

(3.12) 
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(3.14) 


In  order  to  stare  the  iteration  the  state  of  the  structure  at  time  t  is 
used  as  a  first  approximation  to  the  state  of  structure  at  time  tn+^  for  all  the 
terms  except  F^+*.  Since  the  mass  and  stiffness  matrices  change  slowly  with 
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time,  above  assumption  is  adequate  for  the  first  estimation  of  these  terms. 


Because  is  a  nonlinear  function  of  both  gn+*  and  gn+^,  however,  the 


following  predictors  are  used  to  ensure  the  stability  of  solution  [12j. 
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Note  that  Eqs.  (3.15)  and  (3.16)  are  obtained  by  setting  gn+1  =  0  in  Eqs.  (3.6) 
and  (3.7).  For  subsequent  iterations,  the  displacment,  velocity  and  accelera¬ 


tion  are  corrected  by  Eqs.  (3.6)  ~  (3.8). 

n+1 


The  iterations  continue  until  Ajg.  converges  to  within  the  preassigned 
tolerance  A,  that  is, 


a  n+* 

n  aa<  -ii  < 
a r1  -  3" 


where  H*H  is  the  Euclidean  norm  of  '*' 
A  =  0.01. 


For  the  problems  considered  here, 


3.2  Constraint  Conditions 

The  deployment  procedure  adopted  in  the  present  study  models  a  method 
described  in  Reference  11  in  which  a  pushrod  is  connected  to  two  adjacent  bars 
to  transfer  torque  across  the  hinge.  The  deployment  is  accomplished  by  means  of 
driving  a  pushrod  link.  For  this,  it  is  reasonab1e  to  treat  the  relative  angles 
between  two  adjacent  bars  as  a  prescribed  quantity.  Assuming  deployment  in  a 
plane  for  the  sake  of  simplicity,  the  constraint  conditions  are  given  as 

dj  -  qi  =  9  (3.17) 

where 


q-,qj  are  nodal  rotational  angles  with  j  >  i 

6  is  the  prescribed  relative  angle  in  radians 
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Using  Eq.  (3.19),  the  nodal  displacement  vector  3  can  be  written  as 
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Tg:  n  x  n  matrix  with  Tj  ■  =  1,  -  1  and  all  other  components  =0 

3*  =  LqL  ..  qi  ..  9  ..  qnJ 


Since  0  is  a  given  value,  the  virtual  displacement  vector  63*  can  be  transformed 
as  follows; 


'<5q  N 
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53  =  Ib  <  *  )  H  Ib  *3 
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The  equations  of  motion  Eq.  (3.1)  in  virtual  work  form  can  be  written  sym¬ 
bolically  as  follows; 


<5g'  (M  g  +  K  g  -  F)  =  0 


(3.20) 


Substituting  Eqs.  (3.18)  and  (3.19)  into  Eq.  (3.20)  results  in 

53*  (IR  m  Ir  a  +  IR  i<  Ir  g*  ■  Ir  £)  2 


(3.21) 


In  order  to  observe  the  effect  of  deployment  on  the  rotational  motion  of  the 

rigid  mass,  parametric  studies  were  performed.  First,  while  the  total  mass  of 

bars  kept  constant,  the  length  of  bars  was  varied.  Of  course  this  results  in 

changes  in  the  cross-sectional  properties  such  as  area  and  moment  of  inertia 

etc.  Two  different  lengths  were  considered.  They  are  =  300  /2  ft  and 

l_2  =  300  /3  ft.  The  cross-sectional  properties  corresponding  to  L  =  300  /2  are 
“2  2 

area  A  =  4.56  x  10  ft  ,  cross-sectional  area  moment  of  inertia  I  =  Iyy  =  Izz  = 
2.89  x  10”3  ft4.  Next  the  ratio  of  rotary  inertia,  I RAT  10,  was  introduced  as  a 
parameter  as  follows; 

tratto  =  R°tafy  moment  of  inertia  of  fully  deployed  bars 
iK  1U  Rotary  moment  of  inertia  of  rigid  mass 

The  rotary  moment  of  inertia  was  defined  with  respect  to  the  center  of  rigid 
mass.  In  the  present  calculation,  the  values  of  IRATI0  were  chosen  as  1  and  5. 
The  number  of  flexible  bars  were  either  two  or  three. 

Figures  3  to  5  show  the  rotational  angle  of  rigid  mass  vs.  time  up  to  200 
seconds.  In  all  cases,  the  effect  of  flexibility  of  the  bars  are  small  in  the 
initial  stage  of  deployment.  However  the  difference  between  the  rigid  bar  model 
and  the  elastic  bar  model  increases  as  time  progresses.  As  shown  in  the  figures 
the  difference  is  quite  substantial  after  the  completion  of  deployment 
(t  =  100  secs).  For  a  fixed  rotary  inertia  ratio,  the  rigid  mass  with  longer 
bars  (L  =  300  /T)  experiences  a  wider  excursion  of  rotational  angles  than  the 
one  with  shorter  bars  (L  =  300  /IT).  With  a  fixed  L,  the  higher  rotary  inertia 
ratio  ( I RAT  1 0  =  5)  case  shows  more  pronounced  effect  than  the  lower  ratio 
(IRATI0  =  1  case)  as  expected.  For  the  same  I RAT  1 0  and  L,  the  trend  for  the 
three-bar  case  is  similar  to  that  for  the  two-bar  case.  However  the  three-bar 
case  shows  less  effect  of  flexibility. 


4.2  Case  2  -  Deployment  in  Earth  Orbit 

This  problem  is  exactly  the  same  as  that  in  the  previous  section  section 
except  that  now  the  fictitious  space  structure  is  assumed  to  be  in  a  200  mile 
circular  orbit  around  the  earth.  The  folded  bars  are  initially  aligned  in  the 
direction  of  the  velocity  vector  tangent  to  the  circular  orbit.  In  the  calcula¬ 
tion,  the  ratio  of  rotary  inertia  was  fixed  as  I RAT  1 0  =  5.  Instead,  a  mass 

ratio  (MRATIO)  was  introduced  as  follows: 

uoATtn  _  total  mass  of  the  bars 
®  mass  of  the  rigid  body 

The  MRATIO  was  chosen  as  either  0.2  or  0.5. 

The  plots  of  rotational  angle  vs.  time  of  two-bar  case  are  shown  in  Figures 
6  and  7,  which  have  almost  same  trends  as  Fig.  3.  Because  the  gravity  gradient 
tends  to  align  the  structure  with  the  radial  direction,  however,  the  effects  of 
deployment  on  rotational  motion  are  reduced.  It  shows  that,  with  the  fixed 
IRATIO,  the  higher  mass  ratio  (MRATIO  =  0.5)  case  is  subject  to  the  more  gravity 
effect  than  the  lower  ratio  (MRATIO  =  0.2).  Figures  8  and  9  show  the  results  of 
three-bar  case.  They  also  show  the  same  tendency  as  the  previous  case. 

5.  CONCLUSION 

A  finite  element  modeling  and  solution  techniques  was  developed  for  flexible 
space  structures  undergoing  large  timewise  change  in  inertia  and  stiffness  pro¬ 
perties  as  a  result  of  change  in  geometrical  configuration  during  deployment  and 
retraction  maneuvers.  Simple  numerical  tests  and  parametric  studies  demonstrate 
the  validity  of  the  present  formulation  and  the  effect  of  flexibility  on  the 
rotational  motion  of  the  rigid  mass. 

The  present  formulation  is  based  on  the  total  Lagrangian  description  for  the 
description  of  deformation  of  flexible  bars.  In  addition,cross-sections  of  bars 
are  allowed  to  undergo  arbitrary  rotation  in  three-dimensional  space.  This  per¬ 
mits  the  description  of  motion  of  the  flexible  bars  with  respect  to  a  simple 
axis  system  attached  to  the  rigid  mass  or  any  point  in  the  structure. 
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